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In  order  to  provide  diagnostic  data  with  which  to 
identify  and  quantify  the  presence  of  molecular  fragments 
resulting  from  the  destruction  of  highly  energetic  materials, 
we  have  continued  our  ab  initio  calculations  of  excited 
electronic  states  and  radiative  lifetimes  of  the  small 
diatomic  fragments  CO  and  NH.  Such  calculations  are  designed 
to  complement  experimental  spectroscopic  measurements,  pro¬ 
viding  data  which  cannot  be  easily  obtained  in  the  laboratory. 

We  have  identified  particular  excited  states  in  both  the 
CO  and  NH  molecules  which  appeared  to  be  important  as  possi¬ 
ble  origins  of  as  yet  unobserved  band  spectra  or  as  possible 
photodissociation  or  predissociation  channels.  In  addition 
to  characterizing  the  complete  potential  curves  as  a  function 
of  internuclear  separation,  we  also  will  be  calculating 
excited  state  radiative  lifetimes. 

In  CO  we  have  chosen  to  examine  the  four  lowest-lying 
excited  states  in  the  ^T1  and  '*‘Z  +  manifolds,  and  have  devised 
a  series  of  calculations  to  obtain  potential  energy  curves 
and  wavef unctions  at  internuclear  separations  ranging  from 
R=1.5aQ  to  R=20aQ.  The  excited  states  of  ^IT  and  symmetry 

contain  substantial  amounts  of  Rydberg  character  at  short 
internuclear  distances  and  the  configuration  interaction  wave- 
functions  which  we  construct  must  be  capable  of  describing 
valence-Rydberg  mixing. 

First,  a  full-valence  multi-configuration  self-consistent 
field  (FVMCSCF)  calculation  is  performed  for  the  lowest  ^11 
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state.  The  la  and  2a  orbitals  are  always  doubly  occupied  and 
the  valence  orbital  space  consists  of  3o  -  6a,  lit  and  2tt. 

The  MCSCF  orbitals  are  then  used  in  a  single  excitation  config¬ 
uration  interaction  (Cl)  calculation  to  generate  a  set  of 
natural  orbitals,  from  which  one  can  choose  the  lowest  energy 
Rydberg  molecular  orbitals.  The  valence  orbitals  as  well  as 
the  Rydberg  molecular  orbitals  are  then  used  in  a  large  Cl 


calculation  which  includes  configurations  in  the  full  valence 
space,  single  excitations  from  the  full  valence  space  into 
both  the  Rydberg  orbitals  (7o  and  3tt)  and  the  full  set  of 
virtual  orbitals  (8a-32a,  4 tt  —  1 8 tt ,  66),  and  simultaneous  single 
excitations  into  the  Rydberg  and  virtual  spaces.  The  simul¬ 
taneous  single  excitations  are  necessary  to  allow  the  core  to 
relax  to  a  C0+  description  while  the  Rydberg  orbital  is 
occupied.  The  large  Cl  wavef unctions  consist  of  47,036  con¬ 
figuration  state  functions  (CSF's)  for  IT  symmetry  and  26,951 
CSF's  for  ^£  +  symmetry. 

We  calculate  energies  and  wavef unctions  for  4  states  of 
each  symmetry,  and  the  computer  time  (IBM  3033)  required  is  3 
hours  for  the  '*'11  calculation  and  2  hours  for  the  '*T+  calcula¬ 
tion.  In  order  to  stretch  our  computing  dollars  we  carry  out 
these  c.p.u.  intensive  calculations  only  on  weekends  when  the 
rates  for  computer  time  are  25%  of  the  usual  rates.  Because 
of  this  limitation,  progress  is  slow.  However,  from  the  five 
points  which  we  have  computed  it  appears  that  our  calculations 
are  in  excellent  agreement  with  available  spectroscopic  measure- 


merits  which  appear  in  parentheses:  Tg  (A^II-X^E+)  =  8.27eV 

( 8 . 07eV) ,  T  (B1E+-X1E+)  =  10.53eV  (10.78eV),  and  T  (21II-X1E+) 
e  e 

11.27eV  (11.52eV).  The  dissociation  energies  for  the  A^II  and 

xV  states  are  2.77eV  and  10.97eV  compared  to  experimental 

values  of  3.17  and  11.24eV,  respectively. 

The  first  two  states  of  each  symmetry  are  energy  degene- 
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rate  asymptotically,  and  dissociate  to  C(  P)  +  0(  P).  We 
obtain  this  degeneracy  to  better  than  .03eV  at  R=8.0aQ.  The 
third  and  fourth  roots  of  each  symmetry  are  also  energy  dege¬ 
nerate,  separating  to  C('*'D)  +  0(^D)  at  3.23eV  (exact)  above 
the  ground  state  atoms.  Our  excited  atom  asymtotes  are  split 
by  ~.lev  and  lie  -3.6eV  above  the  ground  state  atoms. 

We  have  completed  calculations  of  the  potential  curves 
and  dipole  moments  of  the  four  lowest  states  of  3E~  symmetry 
in  NH  (ground  state,  X3E  ).  There  is  only  one  published  calcu 
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lation  on  the  second  E  state  which  was  not  of  high  quality. 
Our  calculations  indicate  that  the  state  is  repulsive  and  may 
be  an  important  pathway  in  the  photodissociation  of  NH.  The 
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third  and  fourth  E  states  have  never  before  been  character¬ 
ized,  and  appear  to  be  bound  states  which  dissociate  to 
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N (  S)  +  H(2s,2p).  Rydberg-valence  mixing  is  very  important 
for  an  adequate  description  of  these  three  excited  states; 
and  our  atomic  basis  set  of  Slater-type  functions  include 
diffuse  functions  on  both  N  and  H. 

In  much  the  same  manner  as  described  for  CO,  we  obtain 
valence  molecular  orbitals  (2o  -  4o,  Itt)  from  a  full  valence 
MCSCF  calculation  on  the  X^E~,  and  Rydberg  molecular  orbitals 
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from  a  single  excitation  Cl  calculation.  A  number  of  explora¬ 
tory  calculations  were  performed  in  order  to  understand  the 
types  of  conf igurations  necessary  to  obtain  a  balanced  descrip¬ 
tion  of  both  the  ground  and  excited  states,  and  to  allow  the 
states  to  dissociate  properly  asymptotically.  We  concluded 
that  an  accurate  calculation  of  these  four  states  would  have 
to  include  some  double  excitations  for  a  proper  description 
of  the  X  state  bonding,  and  would  have  to  treat  the  Rydberg 
and  valence  orbitals  equivalently  as  much  as  possible. 

In  order  to  include  some  double  excitations  (a  complete 
set  of  double  excitations  would  have  been  prohibitive)  we 
used  the  most  highly  occupied  natural  orbitals  (5a  's,  4  tt's, 
and  2  6  *  s )  obtained  from  a  pairwise  double  excitation  Cl  cal¬ 
culation  for  a  doubles  space.  Our  final  Cl  wavefunction  con¬ 
sisted  of  two  internal  spaces  containing  four  valence  and  four 
Rydberg  orbitals,  and  two  external  spaces  formed  from  the 
doubles  space  and  the  remaining  virtual  orbitals.  The  final 
Cl  wavefunction  has  almost  23,000  CSF’s. 

We  obtain  a  ground  state  energy  which  is  close  (.17eV) 
to  the  full  singles  and  doubles  result  and  a  satisfactory 
De=3.29eV  (3.41eV,  experimental).  Excitation  energies  to  the 
excited  states  are  Tg  =  9.6eV,  10.56eV,  and  11.08eV,  and 
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dipole  moments  from  the  X  T.  state  are  m  good  agreement  with 

those  obtained  from  other  accurate  calculations.  The  asymp- 
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totic  energy  splitting  between  the  N(  D)  +  (H(  S)  limit  and 
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the  ground  state  atoms  N(^S)  +  H(^S)  is  2.64eV  (2.33eV, 

experimental) ,  reflecting  the  fact  that  orbitals  were  opti- 
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mized  for  N(  S)  and  not  N(  D) .  The  asymptotic  energy  split¬ 
ting  for  the  N(^S)  +  H(2s,  2p)  roots  is  10.29eV  (10.20eV, 
experimental) ,  reflecting  the  residual  difference  in  treat¬ 
ment  of  the  valence  and  Rydberg  orbitals  in  the  large  Cl 
wavef unction . 

Publications  and  Talks 

"Oscillator  strengths  for  infrared  vibration-rotation 
transitions  of  CN  and  applications  to  cometary  resonance 
fluorescence",  G.  Lafyatis,  J.  Black  and  K.  Kirby,  to  be 
submitted  to  Astrophysical  Journal. 
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"Potential  Energy  Curves  for  the  Four  Lowest  1  States 
of  NH",  E.M.  Goldfield  and  K.  Kirby,  a  talk  to  be  given  at  the 
fifth  American  Conference  on  Theoretical  Chemistry  in  Jackson 
Hole,  Wyoming,  June  15-June  20,  1984. 


